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O ; Abstract 

^ • CP violation in the system is discussed in the supersymmetric grand unified 

' , theory (GUT) with the right-handed neutrinos. Above the GUT scale, the right-handed 

Q I down-type squarks couple to the right-handed neutrinos. Due to the renormalization 

' group effect, flavor violations in the lepton sector may be transfered to the right-handed 
down-type squark mass matrix, which affects the CP violation in the B decay. Taking 

\ into account this effect, we compare the CP violation in Bd — > V'-^s and B^ — > 

^ • processes. We will find that a significant difference is possible between the CP violating 



phases in two decay processes. 



Many efforts have been made to understand the CP violation. In particular, exper- 
imentally, 5-factories are now being to measure the CP violation in the B system and 
to obtain insights into the physics behind the CP violation. One of the most important 
purposes of the i?-factories is to check the unitarity of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix Vckm using the —>■ ipKs decay mode; in the standard model, the phase 
arg(-[VcKM]cd[V'cKM]c6/[VcKM]fd[VcKM]if') is determined with Bd ipKs Since the phase 
in the CKM matrix is the only source of the CP violation in the standard model,|^ such a 
measurement will provide severe constraints on the CP violation in other processes. 

There can be, however, extra contribution to the CP violation from a new physics beyond 
the standard model. Since the standard model suffers from the so-called hierarchy problem, 
i.e., the stability of the Higgs mass parameter against the radiative corrections, we are 
forced to introduce a new physics at the electroweak scale. Furthermore, the deficits of the 
atmospheric and solar neutrino fluxes strongly suggest non-vanishing masses and mixings in 
the neutrino sector. Since the new physics is, in general, CP violating, the 5-factories may 
be able to observe phenomena related to the new CP violating interactions. 

In this letter, we consider super symmetric grand unified theory (GUT) with the right- 
handed neutrinos as the new physics. In this model, the hierarchy problem is solved by the 
cancellation of the quadratic divergences between bosonic and fermionic loops due to the 
supersymmetry (SUSY), while the neutrino masses are generated by the seesaw mechanism 
0. In such a model, the right-handed down-type squarks interact with the right-handed 
neutrinos and colored Higgs. Of course, such interactions are negligible for the low energy 
physics at the tree level since the colored Higgs is as heavy as the GUT scale. 

However, they become important through the renormalization group (RG) effect 
(Similar effects have been studied for the low energy flavor violating processes 0, |^, ^].) In 
Ref. [^], it was discussed that the right-handed neutrinos affects the structure of the flavor and 
CP violations in the mass matrix of the right-handed down-type squarks in SUSY GUT.p^ In 



particular, SUSY contribution to the e/< parameter can be as large as the currently measured 
value of ex if we adopt an 0(1) phase in the neutrino Yukawa matrix. In this letter, we 
consider different effects, that is, CP violation in the decay processes of the 5-mesons. We 
will see that the SUSY contribution to the decay amplitude of the B^, — >■ <pKs process can 
be large, and that the phase in the decay amplitude can be as large as 0(0. 1) which is 
much larger than the standard model prediction. We will emphasize that phases in the GUT 
Lagrangian play a crucial role for this enhancement. 

Let us first introduce the model we consider. To make our point clearer, we consider a 
SUSY SU (5) model with the right-handed neutrinos. The superpotential of the model is 

Wgut = [Yui, + VP, [Yn],^ + [Y^]^^ + [M^\^ iV„ (1) 

where \E'j, $j, and Ni are 10, 5, and singlet matter fields of SU{b) in i-th generation, while 
'^-'^We neglect the strong CP problem. 

"^^In Ref. 1^, similar effect was also discussed. However, in this article, one of the most important effects, 
i.e., the effect of the phases in the GUT models, was neglected. 
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H and H are Higgs fields which are in 5 and 5 representations, respectively.^ Here, is 
the Majorana mass matrix for the right-handed neutrinos, and for simplicity, we adopt the 
universal structure: 



(2) 



The Yukawa matrix Yu is a complex symmetric matrix while Yd and Yn are complex matrix. 
We take the basis where the Yukawa matrices become 



Yd. Yn = YmVlQl, 



Yu = V^QqYuVq, Yd 
where F's are real diagonal matrices: 

Yu = diag(?/„, yc, yt), Yd = diag(?/d, y,, yb), Yn 
while B's are diagonal phase matrices: 

0Q 



dmg{y 



AQ) 



■AL) -AL) -AL) 

eL = diag(e*^i ,e*^2 \e*<^3 ; 



(3) 



(4) 



(5) 



where phases obey the constraints (p^^^ + (p^^^ + (p^^^ = and + 02^^ + (f)i^^ = 0. 
Furthermore, Vq and Vl are unitary mixing matrices parameterized by three mixing angles 
and one CP violating phase. 

The quarks and leptons in the standard model are embedded in the SU (5) multiplets as 



VI/, ^ {Q, V^Q^qU, QlE},, ~ {D, Q[L}, 



(6) 



where Qi{S, 2)i/6, Ui{3, l)-2/3, A(3, 1)1/3, Li{l, 2)1/2, and £',(1, l)i are quarks and leptons 
in 2-th generation with the SU{3)c x SU{2)l x U{1)y gauge quantum numbers as shown. 
With the embedding (||), the superpotential for the light particles is 



SSM 



Q^ 



UjHu + Qi 



Y, 



D 



Y. 



YnVl 



LjHu 



-M.^NiNi, 



(7) 



with Hu and being the up- and down-type Higgs fields, respectively. In the SU{b) limit. 
Ye = Yd although this relation does not hold for the first and second generations. We expect 
that some mechanism, like higher dimensional operator suppressed by the Planck scale, fixes 
this problem. We assume such a mechanism does not affect the following analysis. In Eq. (|^, 
the unitary matrix Vq becomes the CKM matrix: Vq ~ Vckm- Furthermore, the neutrino 
mass matrix after the seesaw mechanism is given by 



\m 



M, 



sin^ /5 



2M, 



"^^We assume that other chiral muhiplets, in particular, ones which are responsible for the SU(h) breaking, 
do not affect the following argument. 
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where v ~ 246 GeV, and (3 parameterizes the relative size of the Higgs vacuum expectation 
values: tan/3 = {Hy) /{H^). Therefore, with Eq. (H), Vl plays the role of the neutrino 
mixing matrix . In this letter, we mainly consider the neutrino mass matrix suggested by 
the atmospheric neutrino flux deficit p, and also by the large angle MSW solution to the 
solar neutrino problem [0]: 



~ ( 0, 0.004 eV, 0.06 eV 



/ 0.91 -0.30 0.30 \ 
0.42 0.64 -0.64 
0.71 0.70 J 



(9) 



Notice that, due to Eq. (|]), the neutrino Yukawa couplings y^^s increase as M^^ becomes 
larger with the fixed light neutrino masses. With a too large M,^^, the neutrino Yukawa 
coupling becomes non-perturbative below the Planck scale. With Eq. (|^), it happens when 
Mi.^ ^ 10^^ GeV. Therefore, we only consider the cases with M^^ < 10^5 GeV. 

With the embedding (^, all the GUT scale phases drop off from the Yukawa interactions 
among the light fields. However, the phases 0^^^ and (p^^^ remain in the colored Higgs 
vertices. Above the GUT scale Mqut, interaction among N, D, and He is effective in the 
following form 



Wgvt = Ni 



YnVlQl . . D,Hc + --- = J2y,^ [Vl],, e"^f' N,D,Hc + ■■■. (10) 

Although the effect of the colored Higgs is negligible to the B physics at the tree level, 
the interaction ([T0| ) is potentially important since it affects the soft SUSY breaking mass 
parameters through the RG effect. The simplest way to estimate the effect is to assume 
the universality of the scalar mass at the reduced Planck scale M^, ~ 2.4 x 10^^ GeV. 
This corresponds to adopting the model so-called mSUGRA. With this framework, we will 
estimate the possible size of the SUSY contribution to the CP violation in the B system. 

In this approach, the soft SUSY breaking parameters are parameterized by four free 
parameters: the universal scalar mass mo, the universal A-parameter oq, the gaugino mass 
mcs, and the 5-parameter which is determined by the condition for the proper electroweak 
symmetry breaking. With the universal scalar mass at the reduced Planck scale, we run 
all the parameters down to the electroweak scale and evaluate physical quantities using the 
parameters at the electroweak-scale. Of course, the scalar mass parameters may be non- 
universal at the reduced Planck scale, and if so our results can be modified. Since the off- 
diagonal elements of the scalar mass matrix is the most important for our discussion, however, 
we expect that the mSUGRA approach gives us a reasonable and conservative estimation in 
most of the cases. The CP violation in the B system we will discuss will be more enhanced 
if the non- universal contribution is larger than the running effect. If two contributions are 
comparable, the signal can be smaller due to an accidental cancellation. However, such a 
cancellation requires a tuning of the parameters, and we neglect this possibility. 

Next, we estimate the size of the flavor violating off-diagonal elements in the down-type 
squark mass matrix. Relevant part of the soft SUSY breaking terms is given by 



-^^soft = -[mQ]ijQiQ* - [m-^]ijDiDj - {\AD\ijQiDjHd + h.c. 
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-\ (mcsGG + h.c, 



'111 



where Qi and Di are scalar components of the corresponding chiral superfields while G is 
the gluino field. The gluino mass is given by ttlq = \mG3\. With the mSUGRA boundary 
condition, the off-diagonal elements in the mass matrix of D are approximately given by 



. 2 2m ^* 

(3mo + ao)log— 



log 



12) 



GUT 



and hence non-vanishing off-diagonal elements are generated due to the neutrino Yukawa 
matrix. Furthermore, in general, these off diagonal elements have unknown 0(1) phases. 
For a more precise calculation, we solve the RG equation numerically. With uiq ^ and 
oq, the ratios of the off-diagonal to the diagonal elements are approximately given by 



["^|]l3 



6 X 10"^ X 



M, 



1014 GeVy 



2 X 10"^ X 



1014 GeV 



(13) 



where we used the neutrino masses and mixings given in Eq. (^. The neutrino Yukawa 
couplings are proportional to My^, and hence the RG induced off-diagonal elements are 
approximately proportional to M^^. Similarly, the top Yukawa coupling generates those of 
the left-handed down-type squarks: 



;yt [Vckm\u [VcKMlti (3"^o + «o) 3 log — + log 



M, 



GUT 



weak 



(14) 



where M^eak is the electroweak scale. The phases in [Tn'^lij is governed by that in the CKM 
matrix in the mSUGRA case. Numerically, we obtain 



[m|]n 



2 X lO-l 



9 X 10 



-3 



(15) 



Since the off-diagonal elements [m^]j3 



and [^^]i3 are coefficients of AB ^ operators, they 
change the standard model predictions to the mixing and decay of the i?-mesons. In Fig. ^ 
we show the Feynman diagrams contributing to AB = 2 and AB = 1 processes. 

Before discussing the Bd (pKs process, let us consider Bd "ipKs which is the pri- 
mary target of the -B-factories. In the standard model, CP violation in this process is from 
the Bd-Bd mixing (with the standard parametrization of the CKM matrix). Once the CP 
asymmetry in Bd ^Ks is measured, the phase arg(-[V'cKM]cd[VcKM]cfe/[VcKM]td[V^CKM]i;-) 
is determined. 

If there is a new source of the flavor and CP violations, however, this prediction is 
changed. In the SUSY GUT with the right-handed neutrinos, the sdown-sbottom mixing 
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Figure 1: Feynman diagrams contributing to the AB = 2 and AB = 1 processes. The sohd 
hues are quarks with di = d and s, the dashed ones squarks, the wiggled ones with straight 
hues gluinos, and the wiggled line is the gluon. The "dot" • on the squark lines represents 
the mass insertion. For the box diagrams, diagrams with crossing gluino lines also exist. 

through [ml] 13 and [m^]i3 becomes an origin of the B^-Bct mixing. The SUSY contribution 
to the Bd-Bd mixing is discussed in Ref. [|10| where it was pointed out that [m?- ^]i3/[ml 

has to be larger than 0(10^^) to change the standard model prediction significantly. In the 
present case, these ratios are at most 0(10^'^) with the largest possible value of the right- 
handed neutrino mass, M^^ ~ 10^^ GeV. Therefore, we expect that the SUSY contribution 
to the Bd-Bd mixing is smaller than the future experimental sensitivity. Indeed, we calculated 
Ad)^^^ B , the SUSY contribution to the phase in the Bd-Bd mixing amplitude. We found 
that A(f)^^^_^^^ is typically 0(0.1 %) of the standard model contribution, which is too small 
to be seen in the experiments. 

The CP violation in the process Bd 'ipKs is also affected if the phase in the decay 
amplitude is changed. For the process Bd — >■ ipKs, however, the standard model contribution 
to the decay amplitude is at the tree level while the SUSY ones are at the one-loop level. 
As a result, the SUSY contribution to the decay amplitude is negligible. In this scenario, 
the CP violation in the decay Bd — > ipKs is likely to be consistent with the standard model 
prediction. For Bd 'ipKs, the SUSY contribution is relatively insignificant because this 
process has a tree level decay amplitude in the standard model. 

For processes without tree level decay amplitude, the SUSY contribution may be more 
significant. The process Bd — > (pKs is such a process [0, At the quark level, AB = 1 
operators contributing to Bd 4>Ks have a structure like {sb){ss). Such operators are 
induced only at the one-loop level in the standard model. 

The SUSY contribution to the decay amplitude for Bd 4>Ks is from the penguin and 
box diagrams. (See Fig. |l|.) The AB = 1 effective Lagrangian has the following form: 
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+m^.C^T,1s-[7^ tI^^l^'C + {L ^ R) + h.c, 



(16) 



where is the bottom quark mass, G^^ is the gluon field strength, is the SU{3) 



c 



generator and the indices a and b are the color indices. We calculate the SUSY and the 
standard model contributions to the coefficients. In our analysis, we only consider the 
dominant contribution from the squark-gluino loops, and the expressions for the SUSY 
contribution are given in the Appendix. Then, we estimate the decay amplitude: 



M 



{<PK^\C,s\Bd). 



(17) 



In our calculation, no QCD corrections below the electroweak scale are included. Then, we 
adopt the factorization approximation to obtain 



M 



BK 




RR 



^RL " 



1 



DM 



+ (L< 



-C 



RL 



R), 



:i8) 



where is the SU{3)c gauge coupling constant, A^^^; 
used 



3, and the following relations are 



((/)(p^,6^)|sVs10) 

{K%pK)\s''Yb''\B,{pB)) 



F^^'iPB + PkY + F^^'iPB - PkY 



(19) 
(20) 



Furthermore, kdm is an 0(1) coefficient from the hadronization of the chromo-dipole moment 
operator. We estimate kdm using relations derived from the quark model and the heavy 
quark effective theory |jT^: 



2g2 
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- + 0{ml/mi) 



(21) 



where q is the momentum transfer in the gluon line. Typically, 

m 



which gives kdm — 1-2 [12 



|(m| - \ml + m\) 
In our following discussion, we will present results with 
several values of kdm to show the uncertainty related to kdm- 

Now, we discuss the SUSY contribution to the decay amplitude A1^^^^]^o- There are two 
types of contributions, i.e., one proportional to [m^=]23 and the other proportional to [m^ 



Q 



32- 



In the mSUGRA approach, [?7il]23 is generated by the neutrino Yukawa matrix and is pro- 
portional to e^^'^'^i '""^2 The phase in the CKM matrix is independent of tp^^^^s in the SUSY 
SU{5) model. Since the standard model contribution to the decay amplitude A^^^^^^o 
is proportional to [V'cKMltsf^CKMltb, the right-handed down-type squark contribution may 
have, in general, an arbitrary phase relative to the standard model one. On the contrary, in 
mSUGRA, [r?7,^]32 is approximately proportional to 1/t [V'^^Mltsl^CKMltfe- Furthermore, with 
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the model parameters we will use below, |[m^]32| becomes much smaller than |[?«|,]23| since 
|[KjKM]ts| |[^]32|- Therefore, when y^^ ~ y^, the SUSY contribution is dominated by the 
one proportional to [?t?-|^]23 and A^^^^J]^o is approximately proportional to e*^''^3 -^2 ). Of 
course, if we consider different model, ['Ti^]32 may also contribute. 

The CP violation in the decay process Bd (pKs is determined by the sum of the mixing 
and decay phases: 

J total _ /mix _ , j decay 

- +Ml\,, + MTllU. (22) 



The phase in the mixing is universal for the two processes Bd — ipKg and Bd —>■ (pKs- In 
addition, the standard model predicts very small decay phases for these decay modes. As 
a result, in the standard model, (p^B^'^cpKs should be almost the same as the CP violating 
phase measured in Bd ipKs- However, in the model we consider, the decay phase can be 
different in two cases. In order to estimate the SUSY contribution to the decay phase, we 
calculate the quantity 



A0t3.i.o - tan- ( 1 . (23) 



^d-<^i^o - ^'"^^ I , .(SM) 



Notice that l-A^^^^^^ol is almost independent of the phases ^'^'^^'s, and that ^(p'^^'^'^^j^o is the 
maximal possible correction to the decay phase for a given set of model parameters (except 
the GUT phases). Such a maximal value is obtained when the phases are chosen such that 
#)-0(^)^arg[A<|^^^o]+7r/2. 

In Fig. ^ we plot tan[A0^*'^'^^^^^o] as a function of the lightest down-type squark mass 
mjj^. We plot the results with m^j = 500 GeV, oq = and kdm = 0, 0.5, 1, 1.5, and 
2. Let us discuss the behavior of the SUSY contribution. As mentioned in the Appendix, 
there are three classes of contributions: box, color-charge form factor, and chromo-dipole 
contributions. We found that, within our approximation, box and color-charge form factor 
contributions have almost the same size but the opposite sign in most of the parameter space 
we studied. As a result, there is a significant cancellation between these two contributions. 
In particular, with the parameter we used for Fig. ^, an exact cancellation occurs when 
mjj^ ~ 600 GeV. The chromo-dipole contribution is comparable to the others when kdm ~ 1- 
However, this contribution is very sensitive to kdm, and hence the final result strongly 
depends on the value of kdm- From Fig. ||, we see that ^ako can be as large as 0(0. 1) 

with the reasonable value of kdm ~ 1, although more precise calculation of JichK*^ 
requires better understanding of the hadronic matrix elements. 

Here, we should comment on the process /i 67. With the parameters we used, Br(/i 
67) becomes larger than the current upper bound 1.2 x 10'^^ [0 when mj^^ ^ 900 GeV. 
However this is rather a model-dependent statement since all the mixings in Vl affect Br(/i 
67). If the 3-1 element of is ~ 0.1, or if we adopt a non-universal right-handed neutrino 
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Figure 2: tan[A(/)^*'^^^^^o] as a function of the lightest down-type squarks mass with itiq = 
500 GeV, ao = 0, tan/5 = 3, M^^ = 5 x 10^^ GeV, and kdm is 0, 0.5, 1, 1.5, and 2 from 
below. 



mass matrix, Br(/i —>■ 67) may become smaller. In addition, with the neutrino mass matrix 
suggested from the large angle MSW solution with small mass splitting (i.e., so-called "LOW" 
solution) or with those suggested from the small angle MSW or vacuum oscillation solutions 
to the solar neutrino problem (although they are now statistically disfavored 0), Br(/i — > 67) 
is more suppressed. On the contrary, A^^^'^'^^^^o is sensitive only to the 2-3 mixing in the 
neutrino sector. Therefore, we do not exclude the possibility of m^^ ^ 900 GeV in Fig. ^ Of 
course, even with m^^ ^ 900 GeV, A^^'^^'^^^^o can be 0(0. 1). 

Since the standard model predicts almost the same mixing and decay phases in the 
— > ip^s and 4>Ks processes, /^.(p'^g'^^^^Q ~ 0(0. 1) should be an interesting signal. 

Since the uncertainties are expected to be O(10~^) in the standard model calculation [ p!T[| , 
^^%l^(l,K^ ~ 0(0.1) will be regarded as a sign of a new physics beyond the standard model 
if observed. 

So far, we have discussed the decay of the .B^-meson. If the i?s-nieson is available, the 
similar analysis is possible. The discussion is almost parallel to the case of i^^-meson. For 
the process h ccs (resulting in, for e.g., Bg — > ipr], DgDs), a tree amplitude exists and 
the SUSY contribution to the decay phase is negligible. On the contrary, b — >■ sss (which 
induces Bs — > (prj') occurs at the one-loop level and the SUSY contribution can be significant. 
In the standard model, very small CP asymmetries are expected in these processes. Once 
the SUSY contribution is taken into account, however, the decay phase of 0(0. 1) may be 
induced for Bg (prj' since the standard model contribution is one-loop suppressed. 

In summary, the SUSY contribution may change the phase in B^ —>■ <f)Ks- In the future. 
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comparing this process with Bd — >■ "ipKs, we may see a difference in the decay phases in 
these two processes, which cannot be explained by the standard model. This suggests the 
importance to study various decay modes of the 5-mesons. Since the branching ratio of the 
process (pKs is expected to be (9(10~^) [llj, such a study should be challenging at 

the first stage of the present asymmetric B factories. However, at the second stage, or at 
the hadron colliders, more B samples are expected. Then, the study of the decay modes 
with smaller branching ratios is an interesting possibility to look for a signal from the new 
physics beyond the standard model. Therefore, it is desirable to collect a large number of 
the 5-mesons and to study the CP violation in the various decay modes. 



Acknowledgment 

The author would like to thank M. Tanabashi for helpful discussions. This work is supported 
by the Grant-in-Aid for Scientific Research from the Ministry of Education, Science, Sports, 
and Culture of Japan (No. 12047201). 



Appendix: Coefficients 

In this appendix, we present the expressions for the SUSY contribution to the coefficients of 
the Ai? = 1 operators: 

+m,CTT^,r[r, I'WG^ + {L^R)+ h.c. (A.l) 

We only consider the dominant contribution from the squark-gluino loops, and we use the 
mass-eigenstate basis. 

With the soft SUSY breaking terms given in Eq. ([Tl|), the mass matrix of the down- type 
squarks is given bypl 



where /j, is the SUSY invariant Higgs mass (i.e., so-called the /i-parameter) . The above mass 
matrix is diagonalized by the unitary matrix 

PlMpiAB = ml5AB. (A.3) 

Then, the coupling constant for the rfj-d^-gluino vertex is given by 

Xt, = -V2g,[U^]l^e-^^o^ X^^ = -V2g,[Ud]U^e^'^o ^ (A.4) 



"^^In Eq. (A.2) expression, we omit the Yukawa and D-term contributions to the diagonal terms, which 



are included in our numerical calculation. 
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where (pQ is the phase in the gluino mass 

mG3 = |?7iG3|e"^*'^G = m^e" 



(A.5) 



In our calculation, we take = to evade the constraint from the electric dipole moments. 

Neglecting the left-right mixing in the down-type squark mass matrix, the box contribu- 
tion is 



C 



RR 



1 



Box 



EyR* vR vR* vR 
^sA^bA^sB^sB 



Bi{xa,xb) + - 



X 



8N? 



B2{xa,xb) 



c 



(A.6) 



^RL 



1 



Box 



EvR* v-R v-L* vL 
^sA^bA^sB^sB 



X 



iV2 +1 



^ w-Bi{xa,xb) 



AN, 



c 



B2{xa,xb) 



(A.7) 



'^RL 



Box 



~ \ ^ vR* vR v^* 

^g^2^2 Z^^sA^bA^sB^sB 



X 



G AB 

1 „ . . 



ANc 



Bi{xa,xb) + 



ANc 



B2{xa,xb) 



(A.8) 



where the "tilde" denote the SUSY contribution to the coefficients. Here, the functions Bi 
and B2 are given by 



Bi{xa,xb) 



B2{xa,Xb) 



x\ logXA 



x\ logXB 



A{xa - Xb) {xa - 1)^ A{xb - Xj^ {xb - 1)2 

1 



4{XA-1){XB-Iy 
Xa logXA 



Xb log Xb 



(xa- xb){xa-1)^ (xb - xa)(xb - 1)^ 

1 



{xa - 1){xb - 1)' 



(A.9) 



(A.IO) 



(A.ll) 



with Xa = ^j^/'"^!'- Notice that, approximately, the following relation holds in our case: 

Y,X^;X^BBi{xA,XB) ^ Y.^^i^^BBi{xA,XB) ^ 2^|S,(xa, Xg) , 
B B 

where Bi = Bi, B2. 

The contribution from the color-charge form factor are 



C 



RR 



c 



CF 



2Nc 



1 /=^CF r< 



V 
RL 



1 



CF 2Nc 
10 



^R 1 ^RL 



'^R 1 



CF 



(A.12) 



where 



9l 



EvR* -yR 



G A 



(A.13) 



with 



Ci(x) 



2x^ — 9a;^ + 18x — 1 1 — 6 log x 
36(1 ' 
-IQx^ + 45x2 - 36a; + 7 + Qx^{2x - 3) log x 



36(1 -x)^ 



(A.14) 
(A.15) 



For the chromo-dipole operator, 



/=VDM 



93 



E 



(A.16) 



where 



Di{x) 
D^ix) 
Dsix) 
DAx) 



—x^ + 6^2 — 3x — 2 — 6x logx 
6(1 ' 
— + 1 + 2x log X 

(X - 1)3 ' 

2x3 + 3x2 _ g^; _|_ X — 6x2 ^ 
6(1 -x)4 ' 
—3x2 _^ _ X + 2x2 log X 



(A. 17) 
(A. 18) 
(A.19) 
(A.20) 



(x - 1)3 

Other coefficients, Clli C^ji, C'f^, and C^^, are obtained by interchanging L ^ R. 
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